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\ We conduct gravitational microlensing experiments in a galaxy taken from a 

■ cosmological N-body simulation. Hypothetical observers measure the optical depth 

and event rate toward hypothetical LMCs and compare their results with model 
Q ■ predictions. Since we control the accuracy and sophistication of the model, we can 

|H \ determine how good it has to be for statistical errors to dominate over systematic ones. 

22 ' Several thousand independent microlensing experiments are performed. When the 

"best-fit" triaxial model for the mass distribution of the halo is used, the agreement 
between the measured and predicted optical depths is quite good: by and large the 
k>( \ discrepancies are consistent with statistical fluctuations. If, on the other hand, a 

j_j ■ spherical model is used, systematic errors dominate. 

5^ \ Even with our "best-fit" model, there are a few rare experiments where the 

deviation between the measured and predicted optical depths cannot be understood 
in terms of statistical fluctuations. In these experiments there is typically a clump 
of particles crossing the line of sight to the hypothetical LMC. These clumps 
can be either gravitationally bound systems or transient phenomena in a galaxy 
that is still undergoing phase mixing. Substructure of this type, if present in the 
Galactic distribution of Machos, can lead to large systematic errors in the analysis of 
microlensing experiments. 

We also describe how hypothetical WIMP and axion detection experiments might 
be conducted in a simulated N-body galaxy. 
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1. Introduction 

Four years ago the MACHO (Alcock et al. 1993) and EROS (Aubourg et al. 1993) 
collaborations announced candidate gravitational microlensing events toward the Large Magellanic 
Cloud (LMC) demonstrating the viabihty of a new and potentially powerful probe of dark matter 
in the Galactic halo. Microlensing experiments (Paczynski 1986; Griest 1991) are sensitive to 
any object that is smaller than its Einstein radius (in the halo, these objects are known as 
Machos for massive compact halo objects) and therefore complement direct observations, which 
survey the visible content of the Galaxy, and dynamical studies, which measure the total mass 
density. Unfortunately, microlensing experiments are subject to a number of limitations which 
make interpretation of their results rather difficult. First, microlensing events are extremely rare 
and it is necessary to monitor O(IO^) stars in order to get just a few events per year. This has 
restricted present day searches to regions of the sky where there are dense concentrations of stars. 
Indeed the only published microlensing events that can be associated with halo objects have been 
toward the Large and Small Magellanic Clouds (Alcock et al. 1993, 1995, 1996, 1997b; Aubourg 
et al. 1993). Microlensing experiments, with only one or two lines of sight, tell us little about the 
structure of the halo and, more to the point, are sensitive to potentially large systematic errors 
due to our incomplete knowledge of the halo's structure. Moreover, because of the small number 
of events, there are large statistical uncertainties. Finally, present day experiments are unable 
to determine unambiguously the mass and velocity of a given lens. Despite these limitations the 
MACHO collaboration estimates, from the two-year data set, that ~ 50% or more of the mass in 
the halo within 50kpc is composed of ~ 0.5 Mq objects (Alcock et al. 1996). If true, these results 
would have important implications for cosmology, galaxy formation, and star formation. 

The mass density in Machos is determined by comparing the observed number of events 
with the number predicted for a particular model of the Galaxy. Potential systematic errors are 

estimated by seeing how the predicted number of events varies for different "reasonable" Galactic 
models. This strategy, however, is hindered by our limited ability to construct and analyze the 
full complement of acceptable models. In particular, the models used to analyze MACHO's results 
do not generally include triaxiality, velocity space anisotropy, and substructure, all of which are 
expected in the real distribution of Machos. 

In this paper we propose an alternative approach to understanding gravitational microlensing 
experiments. Microlensing experiments are conducted in a halo taken from a cosmological N-body 
simulation. Hypothetical observers make measurements of the optical depth and event rate and 
compare their results with model predictions much as real observers would. We can perform 
a large number of microlensing experiments on a single N-body halo by simply changing the 
positions of the observer and LMC. Since we control the accuracy and sophistication of the model, 
we can determine how good it has to be for statistical errors to dominate over systematic ones. 

Our intention is not to create an N-body realization of the Milky Way Galaxy. Indeed our 
galaxy bears little resemblance to the Milky Way. Moreover, the simulation assumes a cold 
dark matter universe which may not be appropriate if a significant fraction of the dark matter 
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is composed of Machos. However the distribution of particles in the simulated galaxy does 
exhibit the general characteristics (e.g., triaxiality, substructure) that one expects in a realistic 
distribution of Machos. We can therefore test how well theory agrees with observation when these 
features are not fully taken into account. 

Wc find that a simple triaxial model, with the axial ratios and density profile determined from 
the simulation, provides excellent agreement between measured and predicted optical depths. In 
all but a few rare experiments the discrepancies are consistent with the statistical fluctuations that 
one would expect had the particles been chosen at random from the model distribution function. 
If instead an axisymmetric spheroidal model is used, systematic errors become important. And if 
a spherical model is used, systematic errors dominate. 

Microlensing experiments are also sensitive to the velocity space distribution of the Machos. 

If the velocities assumed in the model are too high, for example, one will tend to overestimate 
the mass of individual Machos and their density. To study this, we measure the event rate in the 
simulated galaxy and compare with the expected rate assuming different models for the velocity 
distribution. 

The visible parts of galaxies display substructure such as globular clusters and dwarf galaxies. 
Substructure is also likely to exist in the distribution of dark matter especially in hierarchical 
clustering models of galaxy formation. Microlensing experiments are especially sensitive to 
substructure in the Galaxy: If we are unlucky, a clump of Machos will be passing between us and 
the LMC, biasing our estimates for the density of Machos in the halo (Maoz 1994; Wasserman &; 
Salpetcr 1994; Metcalf & Silk 1996, Zhao 1996). Our simulated microlensing experiment shows 
just such an effect. Even with our best-fit triaxial model, there are a few lines of sight in which 
the measured optical depth is significantly higher than what is predicted. A close inspection of 
these lines of sight reveal that substructure in the halo is often the cause of the discrepancy. 

In Section II we survey the types of models previously considered for the distribution of 
Machos. The methods used to conduct microlensing experiments in an N-body galaxy are 
developed in Section III. Our the results are presented in Section IV. In Section V we describe 
how hypothetical terrestrial dark matter detection experiments can be performed in an N-body 
galaxy. Section VI presents a summary of our results and some concluding thoughts. 

2. Previous Analyses of the MACHO Experiment 

A model for Machos must specify their distribution in mass, configuration space, and velocity 
space. The "standard" halo model, used as a benchmark by the MACHO collaboration, assumes 
that all of the Machos have the same mass, that their velocities are isotropic and Maxwellian, 
and that their density corresponds to that of a cored isothermal sphere (Griest 1991; Alcock et 
al. 1995, 1996). This implies a distribution function (DF) of the form: 



/(x,v) = p(x)F(v) 



(1) 
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where 



A-kG + rl 



(2) 



and 



F(v) 




) 



(3) 



Tc is the core radius, J- is the fraction of the halo in Machos, and Voo is the asymptotic circular 
speed of the total halo. In the standard model, = 5kpc and v^o = 220 km/sec while J- and Mq 
are left as free parameters. It is a likelihood analysis of this 2-parameter model that leads to the 
MACHO collaboration's estimates for the mass and density of Machos in the Galaxy. 

The DF described above represents a highly idealized Macho halo. Deviations between the 
model and the actual distribution of Machos may introduce significant systematic errors in the 
analysis of a microlensing experiment. One can estimate these errors by considering alternative 
models. Several groups, for example, have considered models in which the assumption of spherical 
symmetry is relaxed. For the most part these groups have focused on the seemingly reasonable 
possibility that the halo is an axisymmetric oblate spheroid (Sackett & Gould 1993; Friemann & 
Scoccimarro 1994; Gates, Gyuk, & Turner 1995; Alcock et al. 1995, 1996). Models of this type 
can be constructed by replacing in Eq.(|2|) by m? = + + / where q is the axial ratio 
(< 1 for oblate spheroids) and by multiplying p by \{q) = \/l — q^ / (q'arccosq). Interestingly 
enough, the optical depth and event rate do not change much as q varies from 0.4 (E6 oblate) to 
1.0 (spherical). This is because two competing effects cancel approximately. As the model halo is 
made more oblate, the central density must be increased if v^o is to be kept fixed. This tends to 
increase the optical depth. At the same time our line of sight to the LMC passes through less halo 
material which tends to reduce the optical depth. 

The results for flattened halos have lead to the conjecture that the total inferred mass of 
Machos within 50 kpc is relatively independent of the assumed model ( Gates, Gyuk, & Turner 
1995; Alcock et al. 1996). However there is little direct evidence to suggest that the dark Galactic 
halo is axisymmetric and oblate. Observational clues about the shape of the dark Galactic halo 
come from models of the metal-poor stellar halo (Gilmore, Wyse, &; Kuijken 1989; van der 
Marel 1991), the outer satellites and globular clusters (Hartwick 1996), and HI gas near the 
Galactic plane (Merrifield &: Oiling 1997). The results are highly model dependent and somewhat 
ambiguous leaving open the possibility that the Galactic halo is prolate or triaxial. 

Numerical experiments may provide the best hope for understanding the structure of dark 
halos. The halos found in N-body simulations of dissipationless gravitational collapse are generally 
triaxial with a slight preference to be prolate rather than oblate (Dubinski & Carlberg 1991; 
Warren et al. 1992; see also the simulation discussed below). The inclusion of dissipational matter, 
which tends to settle into a thin disk, will change these results somewhat typically driving systems 
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to become more oblate (Katz and Gunn 1991; Dubinski 1994) though still generally triaxial. And 

there are exceptions (Evrard, Summers, and Davis 1994) suggesting that some halos may be 
prolatc-triaxial even when dissipational matter is included. Witli tliis in mind Holder & Widrow 
(1996) have calculated the optical depth and event rate for prolate halos. For these models the 
central density is lowered relative to what it would be in a spherical model. This effect dominates 
so that predicted optical depth is reduced by a significant amount. 

Galactic microlensing experiments are also sensitive to the velocity space distribution of the 
Machos. If, for example. Machos are preferentially on radial orbits, than the timescale for events 
would be systematically longer than what would be expected assuming an isotropic velocity 
distribution (Evans 1996). The triaxial N-body halos discussed above are supported by anisotropic 
velocity dispersion which can affect the event rate and event duration (Holder k. Widrow 1996; 
Evans 1996). 

One criticism of the models discussed above is that they do not represent true equilibrium 
systems: with the exception of the special case Tc = 0, q = I (the singular isothermal sphere) 
a distribution function given by Eqs. (1-3) does not satisfy the time-independent collisionless 
Boltzmann equation. Evans and Jijina (1994) have attempted to address this concern by using 
the so-called power-law models for the Macho DF. These models are constructed from simple 
power-law functions of the energy and angular momentum and therefore automatically describe 
equilibrium systems. Their main advantage is that they are simple and analytic making lensing 
calculations relatively easy. They do have certain drawbacks. First, while the equipotential 
surfaces are spheroidal, the isodensity surfaces are dimpled at the poles. Second, once the disk is 
included, the models no longer describe self-consistent equilibrium systems. 

Finally, several groups have considered the implications of subgalactic clustering in the 
distribution of Machos (Maoz 1994; Wasserman & Salpeter 1994; Metcalf & Silk 1996; Zhao 1996). 
In general these analyses make ad hoc assumptions about the distribution and structure of Macho 
clusters. 

3. Preliminaries 

3.1. The Simulated Halo 

Our halo is taken from a collisionless N-body simulation of a cold dark matter (CDM) 
universe. Of course Machos are, in all likelihood, baryonic suggesting that a simulation of a 
baryon dominated universe might be more appropriate. However once Machos form they are 
essentially collisionless. We therefore expect that a Macho halo will exhibit triaxiality, velocity 
space anisotropy, and substructure, much like a CDM halo. 

The halo is extracted from cosmological collapse simulation run with a parallel treecode on 
the Pittsburgh Cray T3E (Dubinski 1996). The halo is comprised of roughly 700,000 particles 
within the virial radius though there are only 260,000 particles within the radius of the LMC 
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(Figure |l|). Each particle has a mass Mp = 2.7 x 10^ Mq. The experiments are conducted at 
two different epochs corresponding to ages of 9.7 Gyr and 16.0 Gyr. The results are essentially 
the same for the two time frames and we present only those from the latter. The galaxy is 
prolate-triaxial and has a relatively fiat rotation curve out to large radii (Figure 2). We model 
the mass distribution of the galaxy by assuming that isodensity surfaces are ellipsoidal so that 
p = p{m) where m? = x'^ + u'^/Qi + z'^/Q2 provided we have chosen the axes to coincide with the 
principle axes of the Galaxy. The density profile is fit to an NFW profile (Navarro, Frenk, & 
White 1996): 

p(m) = ^ 4) 

47rgig2 m{m + a^) 



where Mg is equal to the mass inside an ellipsoidal radius m = 3as- Mg, as, qi, and q2 determined 
directly from the simulation. We also use less sophisticated models, specifically, a spheroidal 
model {qi = q2 = 0.5) and a spherical one. The model parameters are summarized in Table 1. 

Ms and depend sensitively on the region of the halo used in the fit. In other words, Eq. (^) 
does not provide a very good global model for this particular galaxy. This is especially true for the 
16 Gyr time frame where the halo is merging with a satellite. Of course, microlensing experiments 
only probe the region of the halo between the observer and the target stars (i.e., between 8.5 kpc 
and 50 kpc) where the density profile is close to a power law. The parameters in Table 1 are chosen 
to fit the density profile in this region. In fact, any number of fitting formulae could be used. For 
example, a cored isothermal sphere (Eq. (|2|)) with Vc = 2.4 kpc, = 1, and Voo = 200 km/sec does 
just as well since the observer is outside the core radius and the LMC is inside the radius where 
the density profile becomes steeper than r~^. In Figure ^ we show the circular rotation speed 
Vc {GM{r)/r)^/'^^ for the halo. Model I, and the "best-fit" cored isothermal sphere. 

A galactic coordinate system is set up in the usual way. Our hypothetical observer is at the 
origin, the center of the galaxy is on the y-axis a distance = 8.5 kpc away, and the z-axis points 
"north". The hypothetical LMC is at galactic coordinates (/, b) = (—32.8°, 281°) a distance L 
from the origin. For simplicity we assume that the hypothetical LMC, in projection as seen from 
the observer, is a circle of radius i?LMC- The foreground of the LMC is a conical volume O with 
base diameter 2iiLMC and height L. 

Different realizations of the experiment are performed by varying the orientation of the halo 



Model 


qi 


Q2 




as (kpc) 


I (triaxial) 


0.55 


0.44 


5.2 


21.4 


II (spheroidal) 


0.5 


0.5 


5.2 


21.2 


III (spherical) 


1.0 


1.0 


6.5 


15.9 



Table 1: Model parameters 
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keeping the observer, LMC, and center of the halo fixed. Since our halo does not have a disk, the 
observer is not constrained to any particular plane. We therefore have three degrees of freedom 
(the Euler angles) in choosing the orientation of the halo. Equivalently, we can imagine keeping 
the halo fixed and considering a family of hypothetical observers positioned on the surface of a 
sphere a distance from the center of the halo. The hypothetical LMCs lie on a concentric sphere 
of radius (L^ — 2rsLcosbcosl + r"^)^^"^ (Figure |l]). For each observer, there is a circle of LMCs a 
distance L away that have the correct orientation. 



3.2. Optical Depth 

The optical depth r is the probability that light from a star in the LMC will be amplified by 
a factor A. This requires that a Macho pass within a distance R = uRe from the line of sight 



where u = 2^^'^ 



1 /2 

A (A"^ — l) — 1 , Re is the Einstein radius: 



o f M L z' (L - z')\'^^^ 
flKM=9.65xlO-(— — ^3^) kpc, (5) 

and z' is the distance to the Macho. For definiteness we set A = 1.34 {u = 1) and L = 50 kpc. If 
all of the Machos have the same mass then r is equivalent to the number of Machos in a tube (the 
so-called lensing tube) whose axis is the line of sight and whose radius is Re- The expectation 
value of T, given an analytic model for /, is 



-1/Ap(x) 



M 



dz'p{z')T:R\{z') 



(6) 
(7) 



where /j(x) = /d^v/(x,v). We also calculate ar, the standard deviation one expects for r 
assuming particles are chosen at random from the model DF. The relevant equations for this are 
derived in the Appendix. 

We can measure r in our N-body halo by counting the number of particles in a lensing tube 
whose radius is given by Eq. (|5|) with M = Mp. (Recall that, for fixed mass density in Machos, 
r is independent of M.) Of course only a very small fraction of the tubes will contain a particle. 
We can mimic what is done in the MACHO experiment by counting the number of particles in a 
large number (~ 10^) of tubes each of which lies, more or less, in the volume However, it is 
more efficient to use all of the particles in weighting them by an appropriate geometric factor. 
To be precise we write the DF for our N-body galaxy as the sum of 5-functions: 



/ = Mp^J(x-x,)5(v-Vi) 

i 



(8) 
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where i labels the particles in the simulation. Substituting into Eq. (^) we have 




where the sum is over all particles in and we use M = Mp in calculating R^. 

An alternative strategy is to replace with a volume 0, whose shape is the same 
as that of a lensing tube. For example, we can choose (7 to be a tube of radius 
i?tubc(^') = 2i?LMC {L — z')/L^)^^^ SO that the maximum radius of the tube is equal to 
the radius of the LMC. In this case all particles contribute equally to r making it easy to visualize 
the phase space distribution of objects in the tube. Furthermore, near the observer, where the 
density of halo objects is presumably highest, is fatter than We therefore expect more 
objects in the tube and hence better statistics. The main disadvantage of this strategy is that 
it distorts the geometry of a realistic microlensing experiment. The angular diameter of Q is 
= 2tan^^ ((L — z')/z') which is far too big near the observer. We will therefore use Q to 
calculate the optical depth. However we use Q to calculate the event rate where it is difficult to 
estimate statistical uncertainties when il. is used (see Appendix and the next section). 



3.3. Event Rate 

In practice the expected number of events, A'^expi is more useful than the optical depth when 
comparing predictions with observations. iVexp is given by 

iVexp = E / --je{t)dt (10) 
JO at 

where t = 2i?E/^ro is the event duration, is the magnitude of the transverse Macho velocity 
in the lensing tube rest frame, e{t) is the detection efficiency, and E is the "exposure" for the 
experiment given in units of "star-yr" . The differential event rate is given by 

dV = — da dvf{x,v)v-uj cos 9 (11) 

where da = REd(pdz' is a surface element on the lensing tube and d'^v = dvz'V-cadv-^dO. Given 
an analytic function for /, the total event rate F is found by integrating over v and a. To 
calculate dT/dt, the event rate as a function of event duration, we change variables from to t, 
differentiate with respect to t, and integrate over the remaining variables. 

If the velocities of the Machos are isotropic and Maxwellian (Eq. (^)) than the expectation 
values for F and dT/dt are 

{T) = ^^^£dz'piz')REiz') (12) 
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and 



dT 
di 



32 



p{z')REiz')exp 



dz' . 



(13) 



An expression for or, the standard deviation for F assuming particles are chosen at random from 
an analytic DF, is derived in the appendix. 



Eq. (11) represents a flux of particles through a surface, something that is difficult to measure 
in an N-body simulation. To circumvent this problem we multiply both sides of Eq. (^) by wdw 
where w is the distance from the axis of the lensing tube: 



wdwdV = -^(i^X(i'^v/(x, v)?;^ cos6Re{z') 



(14) 



Substituting Eq. (|8|) for / and integrating over the volume Vt, yields an expression for F that is 
directly applicable to an N-body galaxy: 



vr 



(15) 



dV/di is estimated by binning the events according to event duration. 



3.4. An Illustrative Example 

We illustrate the techniques developed above in a toy galaxy where the positions and 
velocities of the particles are chosen at random from a known analytic DF. For definiteness we 
assume that this DF is given by Eqs.(l-3) with = 1, Voo = 200 km/sec, and Tc = 5kpc. The 
galaxy is represented by 200,000 Mp = 2.2 x lO^M© mass particles. 1000 individual microlensing 
experiments are performed, each of which assumes a different orientation for the galaxy. To 
be precise, let (a, /?, 7) be the three Euler angles which describe the galaxy's orientation. We 
take 10 even steps in a (from to 27r), cos /3 (from —1 to 1), and 7 (from to 27r). We choose 
Rlmc = 4kpc. This is roughly a factor of two larger than the real LMC, a choice made to increase 
the number of particles in the experimental volume $7. 

From the analytic DF we calculate (r), (F), and {dV/dt) as well as (Jr and crp. We next 
compute, for each experiment, the normalized errors 

= ^ (16) 

er ^ ^ . (17) 

or 
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Figures and ^3 give the probability distributions P{Sr) and -P(er) for the 1000 experiments. 
P{£t) and -P(er) are each well- approximated by a Gaussian of unit variance as they should given 
that the particles are drawn at random from the same distribution used in the model calculations. 

It is also useful to calculate the fractional error 
T- (t) 

= (18) 

liT = . (19) 

/ir and /xr translate directly to the fractional error in the hypothetical observer's estimate for the 
mass fraction in Machos. In this example, where the model is spherically symmetric, (r), (F), 
o"t-, and or are the same in each experiment. We can therefore obtain P{^) from P{e) simply 
by substituting /x^ = t^Orl^r) and ^up = erOT/(F). The situation is more complicated when an 
aspherical model is used since (r), (F), a,-, and or will then depend on the orientation of the 
galaxy. In the analysis below, we use and er to illustrate deviations from Gaussian statistics 
and and [Iy to indicate the potential errors one might encounter in determining the mass 
fraction of Machos. 

Figure ^ shows the measured dV /di for one of the experiments as compared with the predicted 
dV/di from Eq. (111). 



4. Results 

12^ = 1728 microlensing experiments are performed in the N-body halo described in Section 
3.1. The optical depth and event rate are measured and compared with predictions made assuming 
one of the three models of Table 1. The results for our best-fit triaxial model (Model I) are 
presented in Figures |a, |b, and |5|c. Fi gure |5|a is a scatter plot of measured versus predicted 
r for the 1728 experiments. The probability distribution P{£r) for the normalized error is 
reasonably well approximated by a Gaussian of unit variance (Figure ^). There are, however, 
two experiments where the deviations between measured and predicted r differ by more than 4(T^ 
and 14 experiments where the deviations differ by more than Sa-^, significantly more than what is 
expected from Gaussian statistics. The fractional error fir in most of the experiments is < 0.3 
(Figure ^) and comes primarily from statistical fluctuations (cf. Figure ^o). However, in the rare 
experiments where £r ^ 30",-, the fractional error can be as high as 0.7. 

Figures ^ presents the results for the spheroidal model (II). Because of the symmetry in the 
model, different orientations of the galaxy can lead to identical predictions. This accounts for the 
vertical stripes in Figure ^a. Model II is fairly close to the triaxial model: The axial ratios qi and 
q2 have changed by roughly 10%. Nevertheless, there are noticeable systematic effects, as can be 
seen by comparing Figures ^ and |6|c with ||b and |5|c. 
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The results for the spherical model (III) are presented in Figure 0. As in the toy model 
discussed above, the predicted quantities are the same for each of the experiments. We show 
the normalized error distribution P{er)- We can obtain the probability distributions for r 
and fir through the relations r = arSr + (r) and fir = CTrSr/iT) where ar = 6.54 x 10~^ and 
(r) = 5.31 X 10~^. In a large number of experiments the model overestimates the optical depth 
{Sr < 0). These correspond to orientations of the halo where the long axis is near the "galactic 
plane". Likewise, experiments where the model underestimates r find the LMC close to the long 
axis. (This is the orientation considered by Holder &: Widrow (1996).) On the whole, agreement 
between the measured and predicted optical depths is rather poor: The rms normalized error 

= 2.65 corresponding to an rms fractional error of = 0.33. 

Figure ^ presents results for the total event rate F. The triaxial model is used for the mass 
distribution and the velocities are assumed to be isotropic and Maxwellian (Eq. (^)). For this 
model we take Voo = 212 km/sec which is found by computing the rms velocity for all of the 
particles in the halo. 

The experiments with high event rate correspond to orientations of the galaxy which put 
the LMC on the long axis and the observer in the plane containing the two short axes. In these 
experiments, the model overestimates the event rate by as much as 20%. As discussed in Section 
3.1, the parameters in the NFW profile depend sensitively on the region of the galaxy used in 
the fit. The fits in Table 1 are based on the region of the halo inside an ellipsoidal radius of 
50kpc. However, the high F experiments probe smaller ellipsoidal radii where the model tends 
to overestimate the density. A similar problem might arise in actual microlensing experiments 
since the density profile for the Galaxy is determined from observations of stars and gas near the 
Galactic disk while the LMC lies close to the south Galactic pole. 

The triaxial structure of our simulated galaxy is supported by anisotropic velocity dispersion. 
Indeed the velocity dispersion in the x-direction is a factor 1.25 higher than in the y and z 
directions, in agreement with what one expects from the tensor virial theorem for a galaxy of 
this shape. A simple prescription for taking velocity anisotropy into account has the velocity 
distribution in Eq. (|3|) replaced by a modified Maxwellian function of the form (Holder & Widrow 
1996): 



F(v) 



\3/2 



(20) 



^1 and .^2 are related to qi and q2 through the tensor virial theorem. While this model still makes 
the rather dubious assumption that the distribution of Machos in velocity space is independent 
of position it does connect the shape of the halo with its velocity space structure in a reasonable 
way. Unfortunately calculating model predictions with Eq. ( |20|) is quite cumbersome. We choose 
instead to consider a model in which the velocity distribution is assumed to be isotropic but with 

determined separately for each experiment using the particles in the "experimental" volume 
Q,. The results, presented in Figure show a clear improvement over those in Figure ^ where 
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the same Voo was used for all experiments. While ad hoc, this comparison suggests that velocity 
anisotropy can lead to systematic errors in a microlensing experiment. To further illustrate this 
point we measure dT/dt for a single experiment and compare with model predictions. This is 
done in Figure Two models are chosen, one with Voo taken from the simulation as a whole, 
and the other with Voo measured locally (i.e., in the volume fi). Clearly the latter provides better 
agreement between theory and experiment though the net effect on A'cxp (which is calculated by 
integrating over t) may be rather small. 

As noted above, even with the triaxial model, the number of experiments with a relatively 
high > 3cr) discrepancy between measured and observed optical depth is greater than expected 
had the particles been chosen at random from the model DF. We have looked in detail at the phase 
space distribution of particles in these lensing tubes. (Here, we use the volume which has the 
same shape as the lensing tube. This avoids the complication of having particles of different weight 
in the sum for r.) In Figure ^a we show the distribution of particles as a function of v^^ and z', 
the two relevant phase space coordinates for microlensing experiments. A clump of approximately 
50 particles {Mdump ~ 10^ M©), 2/3 of the way to the LMC and moving through the lensing tube 
with a velocity of 450 km/sec, is clearly visible. With enough events, one might "see" such a clump 



in a plot of the event rate as a function of event duration (Figure lib). However, the clump is not 



nearly so pronounced in an "observer's-eye" view of this region of the sky (Figure pT]c). 

Substructure of this type is entirely expected in hierarchical clustering models such as CDM. 



The object shown in Figure 11 is a dwarf halo that will spiral, under the influence of dynamical 
friction, toward the center of the main halo, eventually being stripped apart by tidal interactions. 
Objects of this type arise in about 1% of our experiments, a small but still significant fraction. Of 
course substructure in the stellar halo of the Galaxy has already been observed (see, for example, 
Majewski, Munn, &: Hawley 1994; Ibata, Gilmore, & Irwin 1995). Indeed, observations by Zaritsky 
$ Lin (1997) indicate that there may be an excess of stars between us and the LMC. Though 
these observations are somewhat controversial (Alcock et al. 1997a) they do suggest an alternative 
explaination for MACHO's results: microlensing by foreground stars of a "lumpy" halo (Zaritsky 
& Lin 1997; Zhao 1997). 



5. WIMP and Axion Search Experiments in a Simulated Galaxy 

There are currently a large number of experiments, either in operation or in construction, 
designed to search for elementary particle dark matter candidates such as weakly-interacting 
massive particles (WIMPs) and axions (Dougherty 1995; Rosenberg 1995). Typically these 
experiments measure the density and kinetic energy distribution of dark matter particles passing 
through a terrestrial laboratory. The results of such experiments are therefore subject to 
systematic effects similar to the ones discussed in the previous sections. Clearly, the large scale 
structure of the halo, such as its shape, density profile, and velocity space structure, will determine 
the local density and energy distribution of dark matter particles. Small scale structures such as 
those found in our microlensing experiments can also affect the outcome of dark matter search 
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experiments. Along similar lines, Sikivie, Tkachev, and Wang (1996) have suggested that there 
may be discrete peaks in the local energy distribution of dark matter particles. Their analysis 
is based on the secondary infall model of Fillmore &: Goldreich (1984) and Bertschinger (1985) 
which is characterized by an intricate fine-grained structure in phase space. This model assumes 
spherical symmetry, radial orbits, and smooth accretion of matter and no doubt provides a highly 
idealized description of the Galactic halo. Still the work raises the interesting possibility that 
phase space structure can have an impact on dark matter search experiments. 

We can test some of these possibilities by performing dark matter search experiments in an 
N-body halo. We measure the kinetic energy distribution of particles in a volume V centered on 
a hypothetical observer and compare the results with model predictions. Let be the position 
vector of a hypothetical observer as measured from the center of the halo, p(xs) be the density of 
dark matter in the region of the observer, and -F(v) be the local velocity distribution function, 
i.e., /(xs,v) = p(xs)F{v). The number of particles with kinetic energy per unit mass between k 
and K, + Ak is 



N{k) 



dn , 
— Ak 

an 

p(xjVy2^, 
A 



K,Mp I F{v)d cos 6^c 



(21) 
(22) 



where d^v = v'^dv d cos 6. ud(j)y = y/2KdKd cos 9^ d<j)y. As an illustrative example, we perform 100 
dark matter search experiments taking F to be a sphere of radius 2 kpc. Figure |l2|a shows the 
results for the total number of particles Ny in the volume V . Each point in the plot displays the 
measurement and prediction for a different observer. The model used in making the predictions is 
the same triaxial model discussed above. In addition, we assume that the velocities are isotropic 
and Maxwellian. Figure ^b shows normalized energy spectra as measured by 4 different observers. 
The results are consistent with what one would expect from statistical fluctuations. 

The difficulty we face is that with only 700, 000 particles in the halo, the spatial and energy 
resolution is rather poor. A detailed energy spectrum would require Ny ^ 100 while the structures 
we are interested in may well have a scale significantly less than a kiloparsec. Unfortunately, for 
this simulation, the number of particles in a sphere of radius R is Ny ~ 100 (i?/l kpc)^. 



6. Conclusions and Caveats 



The dark Galactic halo is, in all likelihood, a very complicated system. Most theories of 
galaxy formation would predict that it is triaxial in shape and that the velocity distribution varies 
from point to point and is generally anisotropic. Moreover, if the Galactic halo is built up from 
smaller collapsed objects, as in hierarchical clustering models, it will contain substructure on a 
wide range of subgalactic scales. 
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Traditionally dark matter search experiments such as MACHO have assumed some idealized 
model for the dark matter DF (e.g., Eqs. (1-3)) leading to potentially large systematic errors. 
These errors can be estimated by considering alternative models for the halo such as those 
discussed in Section 2. Our goal has been to test this programme by performing microlensing 
experiments in an N-body realization of a CDM galaxy and comparing the results with predictions 
made assuming a variety of halo models. Our conclusions are as follows: 

• The success of a gravitational microlensing experiment depends crucially on the quality of 
the model chosen for the distribution function of the Machos. The galaxy in our simulated 
experiment is highly asymmetric. By accurately modeling the shape and density profile of 
the galaxy we can achieve excellent agreement between measured and predicted quantities 
such as the optical depth and event rate. If instead, a spherical model is used, the agreement 
is rather poor, the typical systematic errors in the determination of r and F being 30 — 50%. 

• A wide range of fitting formulae will adequately describe the density profile of the halo 
in the region between the observer and the LMC. In this work, we use the NFW profile 
though a cored isothermal sphere, fit to the measured density profile in this region, will 
work just as well. An accurate global fit to the density profile is required only if one plans 
to use observations well outside this region. This suggests that, for actual microlensing 
experiments, one should use a model density profile based on observations of the Galaxy 
between 8.5 and 50kpc. However most of these observations are of material in the Galactic 
plane whereas the LMC is at high Galactic latitude. Systematic errors, such as those seen in 
Figures § and will arise if the halo is far from spherical and the model density profile does 
not provide a good global fit even if we have correctly modeled the shape of the halo. 

• The distribution of Machos in velocity space affects the event rate and event duration both 
of which are used to estimate the mass and density of the Machos. The triaxial shape of 
our N-body galaxy is supported by velocity space anisotropy. By properly taking this into 
account one can improve the agreement between measurements and predictions. 

• Substructure in the distribution of Machos can affect the outcome of a microlensing 
experiment. In particular, a clump of Machos passing between us and the LMC can 
significantly bias the results. Though this situation arises in only ~ 1% of our experiments, 
it can lead to errors as large as 50% in estimates of the optical depth and event rate. 

One should not take the quantitative results in this work too literally. First, the simulation 
was done in the context of a standard cold dark matter universe which may not be applicable if 
Machos are a significant fraction of the dark matter. A baryon dominated universe might be more 
appropriate. Of course, these models tend to have more power on small scales suggesting even 
more halo substructure. Second, our galaxy was created in a pure dissipationless simulation so 
that there is no gaseous or stellar disk. Disk formation will affect the both shape of the halo and 
the density profile. For example, the presence of a disk will probably drive our prolate halo toward 
a more triaxial shape (Katz & Gunn 1991; Dubinski 1994) which might improve the performance 
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of the spherical model. Nevertheless, it is probably optimistic to think that we know the shape of 
the halo well enough to be able to use Gaussian or Poisson statistics for a microlensing experiment. 
Finally, our simulation, with particles of mass Mp = 2.7 x W^Mq, is not able to resolve all of the 
substructure relevant to a microlensing experiment. 

Despite these caveats wc believe our conclusions to be at least qualitatively correct. It will 
be interesting to perform microlensing experiments in other simulated galaxies and in particular, 
simulations that include gas. In a simulated spiral galaxy, for example, one can place hypothetical 
observers in the disk. These observers can then "measure" the rotation curve of the galaxy and 
use this to model the density profile adding another layer of realism to the exercise. Along these 
lines, it would be interesting to choose lines of sight in the simulation toward actual satellite 
galaxies to see if tidal debris of the type discussed by Zhao (1997) and Zaritsky &: Lin (1997) can 
indeed affect microlensing experiments. 
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APPENDIX 

We begin by deriving an expression for the standard deviation cJt-- The expectation value of 
r is given by Eq. (^. In our simulated experiment, we count particles in the conical line-of-sight 
volume Let k be the region in Q such that < z' < + A^.. The probability of finding a 
particle in k is Vk^k = {(^iMCZk)^ ^kp{zk)/]^v- (We assume is small enough so that the 
probability of finding two particles in k is negligible.) A particle in k contributes an amount 
Tk = R^i^k)/ (^LMC-2fc)^ to r. We can therefore write 

{t) =Y.TkVk/^k (1) 

k 

The characteristic function for r is 

^r{s)=X{{l-Vk^k + Vk/^ke'''^'^) . (2) 
k 

Moments of r are found by differentiating with respect to s and setting s = 0. For example 
(r) = ip'^{0)/i in agreement with Eq. (|6|). A straightforward calculation gives 



^ TlVk^k + E E TkVkAkTjVjAj . 

k k jy^k 



(3) 
(4) 



CJt- is given by the expression 



((r-(r)r 



1/2 



dz'T\z')V{z') 



1/2 



(5) 
(6) 



where in the last line we let A^ ^ 0. 

We next calculate the standard deviation err for the total event rate F. For simplicity, we 
assume that the velocity distribution of the Machos is isotropic and Maxwellian. We can write F 
as the following double integral: 

^ ^R^iz')piz')dz' [^^exp(-4-] (7) 
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In this case Vk^k = ''^Rtuhc('^k)^kPizk)/Mp where we are using the volume Cl rather that $7 
(cf. Section 3.2). Similarly the probability that a particle will have v'^ between v^Un and 
Wn + ^n) is V„A„ = exp (— ii„)A„. The characteristic function for T is 

Ms) = nn (l - nAfcV„A„ +PfcAfcV„A„e^G"'=^) (8) 

n k 

where Gkn = "^VoouH^ R^{zk) /ttR^^^^^ The expectation value of F is 

(D = = ^ Vk^kVn^nGkn (9) 

kn 

in agreement with Eq. (j^. The standard deviation of F is 



or 



dz' 



'Re{z')\ p{z') 



\Mtuhe) 



Hubc 
1/2 



1/2 



dz 



Mp 
Mp 



(10) 
(11) 



where Mt^he/MQ = (i?LMc/9.65 x lO^^) (L/50kpc). Had we used $7 rather than 17, we would 
have found a logarithmic divergence. 
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200 kpc 

Fig. 1. — The N-body galaxy in projection. This projection is along the intermediate axis of the 
galaxy so that its prolate nature is readily apparent. The inner circle represents the positions of 
hypothetical observers; the outer circle represents the positions of hypothetical LMCs. 
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Fig. 2. — (a) Density profile and (b) rotation curve for the simulated halo. Crosses represent the 
spherically averaged density and circular rotation speed for the halo. The solid lines are for the 
NFW profile discussed in the text. The dotted lines are for the best-fit cored isothermal sphere 
(rc = 2.4kpc, Voo = 200 km/sec). 
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Fig. 3. — Results of hypothetical mierolensing expermients in the "toy" model galaxy, (a) P{£t) 
as a function of Et', (b) -P(er) as a function of £r- The solid lines are the measured distributions 
for the 1000 experiments. The dotted lines are Gaussians of unit variance. 
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Fig. 4. — Event rate as a function of event duration for a particular microlensing experiment. Tire 
solid line gives the measured r~^dT/di as a function of i. The dotted line is the model prediction. 
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Fig. 5. — (a) Measured versus predicted optical depth r using the "best-fit" triaxial model for the 
galaxy. Each point represents a different microlensing experiment (i.e., different orientation of the 
simulated galaxy), r is given in units of 10~^. There are 1728 experiments in total. In a typical 
experiment, 100 — 200 particles contribute to the optical depth, (b) P{£t) as a function of £r- (c) 
Pi^Hr) as a function of /x^- 
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Fig. 7. — Pi^r) versus Sr for Model III. In this case, the predicted r is the same for ah experiments. 
It is therefore easy to translate this figure into a probability distribution for either r or /Ht- (see 
text). 
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Fig. 8. — Measured versus predicted event rate T. The velocity distribution is assumed to be 
isotropic, Maxwellian, and independent of position. We use Voo = 212 km/sec as determined from 
the halo as a whole. Model I (the triaxial model) is used for the mass distribution, (a) Scatter plot 
of measured versus predicted total event rates for the 1728 experiments, (b) -P(£r) as a function 
of £y- (c) P{ht) as a function of /xr- 
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Fig. 9. — Same as Figure 7 but this time we choose the rms velocity separately for each lensing 
tube. 
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Fig. 10. — Event rate as a function of event duration, dT/di for one of the experiments. The sohd 
Une gives the measured dT/dt. The dotted hne is the model prediction with = 212 km/sec (from 
the velocity dispersion of all of the particles in the halo). The dashed line is the model prediction 
with Voo = 180 km/sec (from the velocity dispersion in this particular lensing tube). 
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Fig. 11. — (a) Phase space distribution of particles where there is a coherent clump in the lensing 
tube. In this experiment the deviation between measured and observed optical depth is over 2>(Tt- 
(a) Distribution of particles in the lensing tube as a function of and z' . (b) dV/di for this 
experiment. The clump would show up as an excess of 70 day events, (c) Map of the sky in the 
region of this hypothetical LMC. The plot gives the angular distribution (measured in radians) of 
all objects in the foreground of the LMC. 
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Fig. 12. — (a) Results for 100 dark matter search experiments. A comparison is made of the total 
number of particles measured in a volume V centered on a given observer with the number predicted 
by the model. The lower and upper solid lines give 2a error bars, (b) Normalized energy spectra 
for 4 of the observers. The dotted lines are model predictions assuming a Maxwellian distribution 
for the velocities. 



